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Abstract 

A complex approach phonon quantum discrete model (PQDM) was developed 
to describe dynamics, kinetics and statistics of phonons in carbon nanostructures 
with zero-chirality of both zig-zag and armchair geometry. The model allows in- 
clude into the pure phonon problem existing interaction with others subsystems: 
electrons, photons, impurities and defects. We predict that planar C- structures are 
geometrically stable and may bridge interelectrode space in strong external electric 
field. The exact solution of generalized thermal conductivity (TC) equation was 
obtained for nanotubes. Temperature distribution along the tube was derived an- 
alytically. The diagonalization procedure for the case ofstrong ph-ph interaction 
was proposed. It was shown the quadratic increasing of heat conductivity with the 
growth of the phonon mean free path (PMFP). Heat capacitance and the entropy of 
carbon linear tubes were calculated as the function of temperature. Our theoretical 
approach explains the nature of good TC in carbon and carbon-like materials by 
existing of the soft vibration branch (low frequency radial breathing mode phonons 
with high density of states at thermal energies) accompanied by structure hard- 
ness (high frequency tp- and z-branches) providing big PMFP. TC coefficient for 
high conducting channel in surrounding medium was calculated. The mechanism of 
heat conductivity temperature damping was analyzed. Two competitive tendencies 
produce TC maximum at intermediate temperatures (100-300)K. It was shown the 
strongly non-linear increasing of effective heat conductivity with the growth nan- 
otubes concentration. It was shown that insertion of armchair nanotube inside a 
medium or its coating by polyacetilene molecule considerably changes the structure 
of radial breathing phonons. 
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1 Motivation 

Phonon dynamics in open and closed carbon and carbon-like nets. How many phonon 
branches are present? Whether exists a simple correlation between differed phonon types 
for different NT geometries? 

Vibrational eigenmode and its amplitudes distribution in the net. 

Vibrational eigenmodes, phonons, micro-sound and heat propagation. Phonon mean 
free path. Sound waves boundary conditions. 

The mechanism of temperature damping at high temperatures and concrete law of 
increasing at low temperatures. 

Maximal thermal conductivity for solitary carbon tube is not measured reliably now 
and changes in wide interval from 200 W/mK to 3000 W/mK for different authors. 

The temperature of maximal thermal conductivity for carbon tubes is not surely mea- 
sured also and may change in wide interval from 150 K to 300 K for different authors. 

The question is whether exist differences in thermal conductivity of solitary carbon 
tubes depending on type (armchair or zig-zag), chirality and diameter. Which physical 
factors are essential in electron transport through the molecule? 

Do phonons (vibrons) participate essentially in electron transport through the molecule? 

Are differences in thermal conductivity of solitary carbon tubes, flat carbon structures 
like graphene and graphite (in plane) essential? 
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2 INTRODUCTION 



Fig.2.1 General picture of connected subjects for the problem of phonon propagation in 
molecular nets. 
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3 Systems under consideration 



Fig.3.1 A polymer molecule absorbed by armchair nanotube. Absorption bonds are shown 
by red lines, C-C bonds of the polymer are shown by violet color. 



Fig. 3. 2 Heat transistor (3-polar contact) on a flat carbon structure connecting three 
baths at different temperatures 7} , T u and T r . 

4 Phonon quantum discrete model (PQDM) 

Jf0 = £^(fcAr + l/2) (1) 



H = £ GUS^S, + bfb iia )+ 

int i h a ^2) 



5 Vibrational dynamics. Linear approximation. 

Zig-zag NT 

mpi = -k'(3pi - pn - p i2 - pa) (3) 
mil = -k(zi - z a ) - 0.25k(2zi - z i2 - z a ) (4) 
mxi = -0.75k(2xi - x iX - x i2 ) - k'(xi - x i3 ) (5) 

Armchair NT 
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mzi = -0.75k(2zi - z a - z i2 ) - k'(zi - z i3 ) (6) 



mxi = -k(xi - x a ) - 0.25/c(2xj - x i2 - x i3 ) (7) 



DA = u 2 A (8) 



u asi A+ s = (C sl , C s2 , C s3 ...C sr ),r = 2n(m + 1) (9) 



Fig. 5.1 The problem's geometry. Z-axes is perpendicular to the figure plane. It's 
shown j and r shifts of an atom. 



Fig. 5. 2. The problem's geometry. F is z-axes projection of the force acting between 
a and b atoms, dz marks the atom a shift along the axes z, dzcos(p/3) is the shift 
projection. 
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Fig. 5. 3. Sketch of the dynamical matrix for zig-zag tubulene created from a graphene 
sheet {5,5} containing 60 atoms by its rolling up. Rose-colored spots are diagonal el- 
ements, green and blue show two types of bonds, rest elements equal to zero, left and 
upper red bands contain atomic numbers. 



Fig. 5. 4. Sketch of the dynamical matrix for armchair tubulene created from a graphene 
sheet {5,5} containing 60 atoms by its rolling up. Rose-colored spots are diagonal ele- 
ments, green and blue show two types of bonds, rest elements equal to zero, left and upper 
red bands contain atomic numbers. 
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6 Frequency spectrum and DOS. Zig-zag geometry. 



g(u sa ) » (10) 

^sa — ^s-l,<7 



Fig. 6.1. Radial mode spectrum of zig-zag tubulene created from a graphene sheet 
{5,15} by rolling up around z-axes. The unity of frequency is ujq 



Fig. 6. 2 Radial mode DOS function for Fig.6.1 system. h=Int[35r/50]. 



Fig. 6. 3 Azimuth mode spectrum of zig-zag tubulene created from a graphene sheet 
{5,15} by rolling up around z-axes. The unity of frequency is u . 
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Fig. 6. 4 Azimuthal mode DOS function (arbitrary units) for graphene sheet {5,15} 
h=Int[35r/50]. 



Fig. 6. 5 Axial mode spectrum of zig-zag tubulene created from a graphene sheet {5,15} 
by rolling up around z-axes. The unity of frequency is u Q . 



Fig. 6. 6 Axial mode DOS function (arbitrary units) for tubulene fragment created from 
a graphene sheet {5,15} by rolling around z-axes. Horizontal axes, frequency in u 
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7 Frequency spectrum and DOS. Armchair geome- 
try. 



Fig. 7.1. Radial mode density of states for armchair tubulene created from a graphene 
sheet {8,15} by rolling up around z-axes. Insertion the same for zig- zag tubulene. 
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Fig. 7. 2 Tangential j-mode density of states for armchair tubulene created from a 
graphene sheet {8,15} by rolling up around z-axes. Insertion is zig-zag equivalent for 
this that is z-mode DOS. h=Int[4r/5]. 



Fig. 7. 3 Axial z-branch density of states for armchair tubulene created from a graphene 
sheet {8,15} by rolling up around z-axes. Insertion is zig-zag equivalent for this that is 
<^-mode DOS. 
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Fig. 7.4 Lower part. Three branches of phonon spectrum for armchair tubulene created 
from a graphene sheet {8,15} by rolling up around z-axes. Upper part. Radial phonon 
band for zig-zag {8,15} NT. 



8 Amplitude distribution and knot's theorem 

Aj = j:\C sj \ 2 n sa , (11) 

S,<7 



Fig. 8.1 Calculated distribution of vibrational amplitudes of {8,8} planar carbon net 
for the 4-th mode. 
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Fig. 8. 2. Calculated r-branch amplitude distribution s=5 along armchair tubulene cre- 
ated from a graphene sheet {6,5} by rolling up around marked direction. Transversal 
knot's lines are degenerated. 



Fig. 8. 3 r-branch amplitude distribution along zig-zag tubulene created from a graphene 
sheet {6,5} by rolling up around marked direction. 5th state. Longitudinal knot's lines 
are degenerated. 



Fig. 8. 4 Calculated mean square amplitude distribution along zig-zag tubulene created 
from a graphene sheet 15,5 by rolling up around marked direction. T=0.03 eV. Circular 
arrows show rolling up of the structure. 

9 Heat transfer 

W^ sa = — IdfgiMN^il + nM) (12) 



| £G i(>s | 2 cdotgi(uj sa ) ■ Ni(u sa ) + \J2 G i r ,s\ 2 ■ g r (uj sa ) ■ N r (u sa ) 
I E G ihS \ 2 ■ gi{uj sa ) + | E G irtS \ 2 ■ g r (uj sa ) 

for 3-polar contact 



n 



s 



(I ^2 G H,s\ 2 ■ 9l{Usa) ■ Ni(Usa) + \J2 G i u ,s\ 2 ' 9u{^sa) • N u (u sa ) + 



+ \^2 G idA 2 ■ 9d{u S a) ■ N d (u sa ))/ 

id 

/(IS Gi ^ s \ 2 ' 9l(Usa) + | Yl G iu,s\ 2 ■ 9u(^sa) + \ G id,s\ 2 ' 9d(^sa)) (14) 



id 
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And for 4-polar contact 



l s = (I J2G kyS \ 2 ■ gi{uj S a) ■ Ni(uj S(T ) + | J2G iryS \ 2 ■ g r {uj sa ) ■ N r (uj sa ) + 

H ir 
+ \J2 G iu,s\ 2 ■ 9u(^sa) ■ N u (LJ sa ) + | J2 G id,s\ 2 ' 9d(^sa) ' N d (u sa ))/ 

/(I Z! G H,s\ 2 ■ 9l(Usa) + I J2 G ir,s\ 2 ■ 9r(u sa ) + \ G iu,s\ 2 ' 9u(^sa) + 

+ \J2 G i d ,s\ 2 -9d(u S a)) 



(15) 



dQ/dt = 2irJ2\ G is\ 2 \ G rs\ 2 u sa gi(uj sa )g r (uj sa ) x 

S,<7 

Ni(uj sa ) - N r {uj so ) 



x 



\ G is\ 2 gii^sa) + \G rs \ 2 g r (uj sa ) 



(16) 



X = L 



dQ/dt 



T r -Tl 



Approximat ion 



- N r (u sa ) » ^^(T r - Ti) 



(17) 



(18) 



\'(L) = 2-kLJ2\ g Is\ 2 \ G rs\ 2 u sa gi(uj sa )g r (uj sa ) x 

ON(uj sa )/dT 



S,<T 



X 



\ G u\giM + \G 

rs\ 9ri^sa) 



(19) 



10 Phonon mean free path and size dependences 



(20) 



C rs — ^ G i r G i, 
The mean free path along the tube lph. 



(21) 



I 



S Id - 1 



ph 



d ' ^ph 

IT' 



In 



t-ph 



<N> 



(22) 



<N>=^ J n{uj)g{uj)duj, 



lph ~ — 



T 



(23) 
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Fig. 10.1. One dimensional phonon dynamics. Calculated total thermal conductivity 
length dependence that includes all vibration branches of armchair NT. Lph is phonon 
mean free path, a is C- bond length. 
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11 Thermal coefficient temperature dependence 

A„ = 2 A (24) 



Fig. 11.1 . Thermal conductivity temperature dependence. Armchair NT. Curve 1 
corresponds to z-branch contribution, curve 2 to j-branch and r-branch contribution is 
presented by curve 3. 



15 



Fig. 11. 2 Thermal conductivity temperature dependence. Zig-zag NT. Curve 1 corre- 
sponds to z-branch contribution, curve 2 to j-branch and r-branch contribution is pre- 
sented by curve 3. 



12 XII. Generalized phonon transport equation 



W^ sa = — \G ls \ 2 g^MMil (25) 



^ = ^ (26) 



Df = \G,„\ 2 gM (27) 



Gisa — G rsa — ^Gi ; Cj iS(7 (28) 

k 

Temperature distribution along the tube 



T(z) = u sa /Ln 1 + — — — r , — (29) 

\ ni{uj sa ) + [n r (uj sa ) - ni{(jj sa )\z J 
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Fig. 12.1 Phonon jumps over a carbon armchair structure (arrows). Nanotubes are 
ideal one-dimensional systems as to the phonon dynamics. 



Fig. 12. 2. One dimensional phonon dynamics. Exact temperature distribution along 
the tube given by (4.4). Tl ,Tr.are temperatures of baths connecting with tube ends, L is 
the tube length. End points come together all modal (partial) temperatures for each ss. 



17 



13 Summary and discussion 



Fig. 13.1. Averaging procedures and speed of propagation, a is interatomic distance, \xi 
is mean atomic shift , x is extended coordinate of the freedom's degree, V is excitation 
velocity, m is atom mass, k is elasticity coefficient. 



Fig. 13. 2. Temperature dependence of thermal conductivity coefficient. Two competi- 
tive tendencies. 

The boundary conditions for microsound spreading through periodic NT complexes 



A complex approach PQDM was applied to describe phonon dynamics and heat 
spreading processes in carbon NT-solid medium composite material and in nanotubes 
embraced by a polymer molecule. 

Atomic vibration dynamics was considered for carbon zero-hirality nanotubes con- 
nected with solid matrix. Vibrational eigenmodes, density of states and amplitude distri- 
bution for tubes inserted into a cylindrical well inside the medium and bounded with its 
hard inner wall were calculated. Phonon band parameters were analyzed in linear approx- 
imation for three types of vibration: athimuthal or tangential j-branch, radial breathing 
mode (r-branch) and longitudinal z-branch. 

It's shown that phonon propagation in actual nanotubes is characterized by a kind 
of " compactification" of circular freedom's degree due to the big phonon mean free path. 
Nanotubes of actual diameters are ideal one-dimensional phonon qnd heat conductors. 




mui 2 < z >= const 



(30) 



< z >= Ae ifz + Be~ lfz 



(31) 
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Phonon band structure was investigated for armchair nanotubes on the base of hier- 
archical law and system symmetry. 

Thermal fluxes and thermal conductivity were considered in PQDM. Temperature de- 
pendences were obtained. The mechanism of heat conductivity high temperature damping 
is cleared. Two competitive tendencies produce thermal conductivity maximum at inter- 
mediate temperatures (100-300)K. 

The exact solution of generalized thermal conductivity equation was obtained for nan- 
otubes. Temperature distribution along the tube was derived analytically. 

Size dependences were considered for thermal conductivity. It was shown the linear 
increasing of heat conductivity with the growth of the phonon mean free path. 

The explanation of the nature of good thermal conductivity in carbon and carbon- 
like materials by existing of the soft vibration branch (low frequency RBM phonons with 
high DOS at thermal energies) accompanied by structure hardness (high frequency j- and 
z-branches) providing big mean free path for phonons (Fig.5.2). 

Adding of new layers or new walls to single-walled NT makes breathing r-branch of 
vibrations harder and causes the sharp decreasing of phonon density of states at the same 
phonon mean free path 

Pressure decreases thermal conductivity. The effect is connected with total hardening 
of all bonds and phonon modes going away from active thermal zone. 

Melting decreases thermal conductivity by another reason: the phonon's mean free 
path becomes small. 

Modified Fourier equation of heat propagation in microscopic molecular channel was 
obtained. 

Boundary problem for temperature distribution in channel-medium system was con- 
sidered in cylindrical symmetry. The method of partial solution building for cylindrical 
boundary problem was proposed. Temperature distribution isotherms were calculated. 

Concentration dependences were considered for thermal conductivity of NT-solid ma- 
trix composite. It was shown the strongly non-linear increasing of effective heat conduc- 
tivity with the growth nanotubes concentration. 

Phonon dynamics was considered for carbon zero-chirality nanotubes embraced by 
polymer molecule. Vibrational eigenmodes, density of states and amplitude distribution 
for a bounded tube were calculated. Phonon band parameters were obtained in linear 
approximation for three types of vibration. 

It was shown that insertion of armchair nanotube inside a medium or its coating by 
polyacetilene molecule considerably changes the structure of radial breathing phonons 
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